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Abstract 

In the present article we analyze the problem of a relativistic Dirac electron 
in the presence of a combination of a Coulomb field, a 1/r scalar potential as 
well as a Dirac magnetic monopole and an Aharonov-Bohm potential. Using 
the algebraic method of separation of variables, the Dirac equation expressed 
in the local rotating diagonal gauge is completely separated in spherical coor- 
dinates, and exact solutions are obtained. We compute the energy spectrum 
and analyze how it depends on the intensity of the Aharonov-Bohm and the 
magnetic monopole strengths. 
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I. INTRODUCTION 



The Dirac equation is a system of four coupled partial differential which describes the 
relativistic electron and other spin 1/2 particles. Despite the remarkable effort made dur- 
ing the last decades in order to find exact solutions for the relativistic Dirac electron the 
amount of solvable configurations is relatively scarce, being the Coulomb problem perhaps 
the most representative example and also one of the most discussed and analyzed problems 
in relativistic quantum mechanics. Among the different approaches available in the litera- 
ture for discussing the Dirac-Coulomb problem in the presence of other interactions like the 
Aharonov-Bohm field or any other electromagnetic potential we have the quaternionic ap- 
proach proposed by HautotB, the Stackel separation method developed by Bagrov et aB, the 
algebraic method of separation of variablesi'i, the shift operator method!, and the algebraic 
method proposed by Komarov and Romanovai. 

Recently, Lee Van Hoang et afl have solved the Dirac-Coulomb problem when an 
Aharanov-Bohm and a Dirac magnetic monopole fields are present. The authors use, for 
tackling the problem, a two dimensional complex space which results after applying the 
Kustaanheimo-Stiefel transformationi on the three space variables, reducing in this way the 
Kepler problem to an oscillator problem. This idea lies on the utilization of a SU(2) dynam- 
ical algebra for computing the resulting energy spectrum!!, which, like the spinor solution, 
is expressed in terms of intrinsic coordinates appearing after using the complex space. The 
utilization of different techniques for studying the Dirac-Coulomb field in the presence of an 
Aharonov-Bohm field or a magnetic monopole could give rise to the idea that this problem 
is not soluble without introducing new variables or additional conserved quantities. Here it 
is shown that using the algebraic method of separation of variables it is possible to solve 
the Dirac equation in the presence of a Coulomb field and a scalar 1/r potential with an 
Aharonov-Bohm and a magnetic monopole fields. The advantage of this approach is that 
does not require the introduction of non bijective quadratic transformations, also it becomes 
clear the role played by the Dirac magnetic monopole as well as the Aharonov-Bohm field 
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in the angular dependence of the spinor \I/(r) solution of the Dirac equation. 

The article is structured as follows: In Sec. II, applying a pairwise scheme of separation, 
we separate variables in the Dirac equation expressed in the local rotating frame, we separate 
the radial dependence from the angular one. In Sec. III. the angular dependence is solved 
in terms of Jacobi Polynomials. In Sec IV, the separated radial equation is solved and the 
energy spectrum is computed. In Sec. V, we discuss the influence of the Aharonov-Bohm 
field and the magnetic monopole on the energy spectrum. 



In this section we proceed to separate variables in the Dirac equation when a Coulomb 
field, a scalar 1/r potential as well as a Dirac magnetic monopole and a Aharonov-Bohm field 
are present. For this purpose, we write in spherical coordinates the covariant generalization 
of the Dirac equation 



where 7 M are the curved gamma matrices satisfying the relation, {7^,7^}+ = 2g M1/ , and 



II. SEPARATION OF VARIABLES 




(1) 



are 



the spin connections^, with V(r) as the scalar 1/r field 




a' 



(2) 



r 



where a' 



is a constant, and the vector potential reads 



A - /t( mon ) _i_ 4(Coul) , 4(AB) 



(3) 



where the components of the Coulomb potential A^ oul ^ take the form 




(4) 



r 



the Aharonov Bohm potential A| t AB - ) reads 



and the Dirac monopole field A^ 10 ^ is 

A(m , = (i^), (6) 

r sin if 

where, following the Dirac0 prescription for quantizing the magnetic charge, g takes integer 
or half integer values 

g= 3 -,j = 0, ±1, ±2.... (7) 

If we choose to work in the fixed Cartesian gauge, the spinor connections are zero and 
the 7 matrices take the form 



7° = 7° = 7°, 7 1 



(7 1 cos ip + 7 2 sin ip) sin if + 7 s cos 



f = - 
r 



1 r^,i 



-2 



(7 1 cos ip + 7 2 sin y>) cos if — 7 s sin if = — , (8) 



7 3 = — : — -(— 7 1 simp + 7 2 cos<^) ~ 



r 

3 



r sin ?9 r sin ^ 

where 7° are the standard Minkowski gamma matrices, and the Dirac equation in the fixed 
Cartesian tetrad frame (§) takes the form 

{f(d t + iV{r))+fd r + ^ + ^^(d v -iF-i{l- C osd)g)+M+V{r)\^ Cart = 
[ r r sin if J 

(9) 

where we have introduced the spinor ^cart, solution of the Dirac equation in the fixed 
tetrad gauge. In order to separate variables in the Dirac equation, we are going to work in 
the diagonal tetrad gauge where the gamma matrices 72 take the form 

7 d ° = 7°, f d = 7\ ll = -7 2 , ll = ^7 3 (10) 

Since the curvilinear matrices 7^ and 7^ satisfy the same anticommutation relations, they 
are related by a similarity transformation, unique up to a factor. In the present case we 
choose this factor in order to eliminate the spin connections in the resulting Dirac equation. 
The transformation S can be written as0 



s 



r(sin^) 1 / 



T7^ ex P(-^7 1 7 2 )exp(-^7 3 7 1 )a = S a 



(11) 



where a is a constant non singular matrix given by, a = |(7 1 7 2 — 7 X 7 3 + 7 2 7 3 + I) which 
applied on the gamma's acts as follows 



37^ 1 = 7 3 , a 7 2 a 1 = 7 1 , a 7 3 a 1 = 7 2 , 



(12) 



the transformation S acts on the curvilinear 7 matrices, reducing them to the rotating 
diagonal gauge as follows 



ST^S = gf'f = 7^ (no summation) 



(13) 



then, the Dirac equation in spherical coordinates, with the radial potential V(r), in the local 
rotating frame reads 

\j°(d t + iV{r)) + 7 1 <9 r . + +—d* + -^s(d v -iF- z(l - cos%) + M + V(r)\ V rot = 
[ r rsmv J 

(14) 



where we have introduced the spinor ^f rot , related to ^cart by the expression 



^Cart = S^ rot = S a^ 



rot 



(15) 



and 7^ are the standard Dirac flat matrices. 



Applying the algebraic method of separation of variables™, it is possible to write eq. (14) 
as a sum of two first order linear differential operators Ki , K 2 satisfying the relation 



K 1 ,K 2 ]=0, {k 1 + K 2 }<5> = 



(16) 



K 2 <& = k$ = -K-lQ 

then, if we separate the time and radial dependence from the angular one, we obtain 

„3 



(17) 



7 



sin v 



7°7 1 $ = /c$ 



(18) 



-ir 



7°<9 t + 7 1 9 r + M + ij°V(r) + V(r)} 7° 7 1 $ 



(19) 



with 



rot 



(20) 



where the operator Ki as well as K2 have been chosen to be Hermitean, therefore the 
constant of separation k appearing in and (|1£]) is real. Notice that if we drop out the 
Aharonov-Bohm and the magnetic charge contributions in eq. ( |T8D we obtain the Brill and 
Wheeled angular momentum K, and separation constant k takes integer values. 

Now we proceed to decouple the equation (|Tj|) governing the angular dependence of 
the Dirac spinor In order to simplify the resulting equations, we choose to work with the 
auxiliary spinor $ related to $ as follows 



(21) 



Since the operators K2 and K\ commute with the projection of the angular momentum 
—id,/,, with eigenvalues m, we have that equation (|1^) takes the form 



7 



siwd 



[m + F + (1 -cosi%*) 



7 ° 7 3 $ 



-ik$ 



(22) 



In order to reduce the equation (2~2~) to a system of ordinary differential equations, we 



choose to work in the following representation for the gamma matrices,^ 



r 



(0^ 


,7° = 


(-i 


"I 












{ 


V 



(23) 



Then, substituting (|2~3"D into (|^) we obtain, 



<r 2 d$ — i 



a 



sin?? 



a 



[m + F + (1 -cos$)g) 



$1 = -ik$i 



& d# + i- — -(m + F + (1 - cos$)g) 



sin $ 



-ik§-) 



(24) 
(25) 



with 



$ = e 



imcj) 



(26) 
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III. SOLUTION OF THE ANGULAR EQUATIONS 



In this section we are going to solve the systems of equations (H) and (^)- It is not 
difficult to see that the spinor $2 can be written as G(r, t)a 3 $i, where G(r, t) is an arbitrary 
function, this property allows us to consider only the system (^). Using the standard Pauli 
matrices , we have that (f24|) reduces to 



ksiwd^i — (m — F — (1 — cos-^g)^ + sini? 



d-d 



— fcsini9£2 + (m — F — (1 — cosi?))£i + sin^—— = 

ait 







(27) 



(28) 



with 



$1 = Qe 



im<j> 



V 6 / 



(29) 



where Q is a function depending on the variables t, r, to be determined after a complete 
separation of variables. In order to solve the coupled system of equations (|2"7|)-(|2"8"|) we make 
the following ansatz, 



£1 = (sin^(cos^)V(^), 6 = (sin^) c (cos^g(tf) 



(30) 



where a, b, c, and d are constants to be fixed in order to obtain solutions of the governing 
equations for /(■#) and g($) in terms of orthogonal special functions. Substituting (|30| ) into 
(f27|) and (|28"D we obtain, 



7 / \ i m F a 1\ /■/ \ i-i \df(x) 

- kq(x) - ( 1 h -)f(x) + (1 - x)^—^ 

HK J K 2 2 2 2 JJK J v 1 dx 



i r, \ ,d m F ,dq(x) 

kf(x) + (--g--- -)q(x) + (1 + x)-±-> 







(31) 



(32) 



where we have made the change of variable x = cos and we have simplified the resulting 



equations (|3T| ) and (|32| ) by imposing c = m — F, and 6 = m — F — 2g. If we set a — m— F + l, 
and d = m — F + l— 2g,and make the change of variables u — (1 — x)/2, we obtain that the 



coupled system of equations (|3T[) - (p2|) takes the form 



kq{u) _ ( m _ F + _ = o (33) 



fc/(«) + (m - F - 2g + l -)q{u) + (u - 1)^ = (34) 



from (13~3T) and (Bl) we obtain 



,„ ,d 2 q(u) ^ 1. _ 1. ,dq(u) 

<l - u )^r + [{m-F + -)-2(m-F-g+ ^>}^ L + 

+ [k 2 -{m-F+^)(m-F-2g + ^)}q{u) = (35) 

s d 2 f(u) „ 3. . „ . ,df(u) 

U{1 ~ U) ^y + [(m-F + -)-2(m-F-g + 1)«]^+ 

+ [fc 2 - (m - F + - F - 2g + ^)]/(«) = (36) 

the solution of the equation Q3"5D can be expressed in terms of the hypergeometric function 
F(a,b; c- u)M as follows 

q(u) = C\F(a, b; c; u) (37) 
where ci is a constant, and a,b, and c are 

a = m- F - 9 + ^~ \/k 2 + g 2 (38) 



b = m-F-g + - + ^/k 2 + g 2 (39) 



c = m -F + - (40) 
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then, with the help of (|34"D and ( [37]) we find that f(u) reads 



/(«) = ci — lrrr^(o> & ; c + 1; M ) ( 41 ) 

Since we are looking for normalizable solutions according to the product 

r 2n + 

2tt / d$ = 5 kk , (42) 



we have that the series associated with the hypergeometric function F(a, b; c; u) should be 
truncated reducing it to polynomials. This one is possible if a = —n or b = —n in ( p?|) where 
n is a no negative integer value. Then using the relation between the Jacobi Polynomials 
and the functions F(a, b; c; x) 

Pt P \x) = r( " + a + 1) F(-n, n + a + P + l;a + l; ^) (43) 
n\T(a + l) 2 



we have that fl4T|) and (^) reduce respectively to 

fix) = c 'p(m-F+l/2,m-F-2 g -l/2) ^ (44) 

= F-l/2,m- F-2 S+ l/2) (x) (45) 

where d is a constant, and n reads 



n= -m+F+g- - + ^k 2 + g 2 (46) 
Then, the components £i and £2 of the spinor $ (]29f) can be written as 

6 = C '(sin ^)™"^+ 1 (cOS f-2 S p j (m-F+l/2 ) m-F-2 S -l/2) (c0S1?) (4?) 

£ 2 = C ^ + v/ f + ^ (sin -) m - F (cOS ^)--^-29+lp j (m-F-l/2 ) m-F-2 S+ l/2) (c0S1?) (4g) 

rC 2 2 

Notice that the orthogonality relation for the Jacobi Polynomials 

1 (1 - x) a (l + xfP^\x)P^\x)dx = 
-1 

2«+^+i T(n + a + l)T(n + (3 + 1) ^ 



2n + a + /? + 1 T(n + l)r(n + a + /3 + 1) 
imposes some restrictions on the values of m, F and g in (f47|) and (|4*8|) In fact, from (|49[) 
we have that a > —1, /3 and consequently we have that 

m- F + ^> 0, m-F-2g + l/2>0 (50) 

affords the required condition of orthogonality. However, some values considered in the 
inequalities given by (|50|) fail in fulfilling a condition which should be also took into account 
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in selecting the possible values of m, F and g. : All the expectation values associated with 
the separating operators should exist. In fact, if we consider the operator K 2 defined by the 
relation (|T8| ) we have that the integrals / <M£i9#£i as well as / dd^d-d^ should be finite. 
This one imposes some restrictions on m, F and g. Looking at the convergency when d — > 0, 
and i? — > it we obtain the relations 

m-F>0, m- F - 2g > (51) 



which are weaker than those given in (|50| ) 

In order to be able to consider other relations among the parameters m, F and g different 
from (fUD we are going to consider a second solution of the equation differential equation 



q(u) = u (1 - M) c - a - fe F(l - a, 1 - 6; 2 - c;tt) (52) 

then, using the recurrence relations for the hypergeometric functions, we find that a solution 
for the system (|33|)-(|34|) reads 



f(u) = c 2 u~ 1/2 - m+F (l - u)- m +f+*9+V* F (i - a, 1 - b; - - m + F; u) (53) 

q( u ) = c 21 - u -^ 2 ~ m+F (l - u ym+F + 2 a+ i/2 F{l _ x _ b 3 _ ^ + F _ 

2 — m + F 2 



where c 2 is a constant. From (p3|) and ( p4|) we obtain that £i and £ 2 take the form 

& = C(sin ^"^(COS !)-^+%+lp(-V2-™+F,2 9 -m+F+l/2) (cos1 5 ) (55) 

£ 2 = c k (sin -) 1 -™+^(coS t)-™+F+2 9 p(l/2-rn + F,2 9 - m+ F-l/2),^ (5g) 

\fh 2 + g 2 — g 2 2 
with c as a constant of normalization, and in the present case n is 



n = m - F- g-- + ^k 2 + g 2 (57) 

Following the same reasoning used for deriving (|5"T|), we have that the expressions for £i and 
£ 2 given by fl55|) and ([56D are valid when 
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-m + F > 0, 2g - m + F > (58) 
A third possible solution of the equation fl35|) can be written as 

q(u) = (1 - u) c - a - 6 F(c - a, c - 6; c; u) (59) 
then, substituting ( |59|) into fl34|) we find that f{u) and g(w) take the form: 



f( U ) = C 3 (l - u )-m+F+2 9+ l/2 F ^ +g+ ^ k 2 +g 2 jl+g _^ k 2 +g 2. m _ F+ _. u) (60 ) 

q(u) = -c 3 m ' F k +1/2 (l - n)-" 1 ^- 1 / 2 ^^ + <? - m - F + 1; u) 

(61) 

where C3 is a constant. Using the relation between the hypergeometric function F(a, b; c; x) 
and the Jacobi Polynomials P^ a,/3 ^(x) given by eq. (H3I) we arrive at 



6 = C(sin ^—^(COS ^-m+F+2 9 +^m-F + l/2,- mW 2 9+ l/2) (cos ^ (62) 

= -c^(sin |) m - F (cos ^-^^p^-iA-^^-i/^^^^ (63) 
where c is a constant of normalization and n is given by 



n = ^k 2 + g 2 -g (64) 
In this case, we have to impose the following restrictions on the values of m, F and g. 

m - F > 0, -m + F + 2g>0 (65) 
Finally, considering as solution of the equation fl35|) the expression 

q{u) = u l ~ c F(a - c + 1, b - c + 1; 2 - c; u) (66) 
then, in the present case the functions f{u) and q{u) read 



/(„) = u- l ' 2 - m+F F{-g - + 0*, -g + ^k 2 + g 2 ;--m + F; u) (67) 
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</< " i = i _ k - ^ 2 - m ^ F F{-g - y/k? + g 2 + 1, -g + v/fc 2 + 9 2 + 1; ~ - m + F; w) (68) 



& = C4 (sin ^)- m+F (cOS l r -F~2 9p (-l/2- m+ FA/2 +m ~F~2 g ) icos ^ ( g g) 
& = C4 ^( sin +^( C0S ^ )m -F-2 9+ l p (l_/2- m+ F,-l/2 +m -F-2 9 ) (cos ^ (7Q) 



n = Jk* + g* + g (71) 



with n given by 



and C4 is a constant of normalization. The solutions (p9|) and (|70[) are well behaved according 
to the Dirac inner product as well as to the expectation value of the operator of angular 
momentum (1181) if 



-m + F > 0, m - F - 2g > (72) 
then, we have that the results flB3j ) and ([71]) can be gathered as follows 



n=jk*+g*-\g\ (73) 
when the condition on m, F and g 

F + g+ | g \> m > g- \ g \ +F (74) 

is satisfied. 

Regarding the eigenvalues m of the projection of the angular momentum operator —id v 
we have that since the transformation (|ITD , relating the Dirac spinors ^ rot and ^cart m the 
local (rotating) and the Cartesian tetrad frames transforms after a rotation as follows 

S t (<p + 2Tr) = -S g (<p) (75) 

and the spinor ^cart is single valued, then we obtain 

^ ro t{ip + 2tt) = -^rot(^) (76) 

and therefore m takes half integer values 

m = N + -, N = 0, ±1, ±2... (77) 
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IV. SOLUTION OF THE RADIAL EQUATION 



Now, we are going to solve the system of equation (|19D , governing the radial dependence 
of the spinor ty rot solution of the Dirac equation. This equation can be written in the form, 



-7 3 <9t + 7°<9r + (M + K(r)) 7 °7 3 - ij 3 V(r) + i- $ = 

r / 



(78) 



Using the representation for the gamma matrices given by (P3"|), and the fact that eq. 
commutes with the energy operator id t with eigenvalues E, we obtain the following system 
of equations 



d r + -)$!- a 3 (E - V(r) - M - V(r))$ 2 = 



(79) 



-dr + - $ 2 - a 3 (E - V(r) + M + V(r))$i = 
r / 



(80) 



From ([24]) - (f25p , (fT9|) - flSfp and the fact that the Dirac equation (|14]) commutes with — ic^, and 
i(9 t , we have that the spinor l> can be written as follows 



$ = c e 



i(rrup—Et) 



1 £iA(r) ^ 
6A(r) 
cCiB(r) 
-c£ 2 B(r) 

where c is a constant, and A(r) and B(r) satisfy the system of equations 



V 



/ 



d r + - ) A(r) - (£7 - V(r) - VYr) - M)B(r) = 
r / 



(81) 



(82) 



d r + - J B(r) - (£? - V(r) + V"(r) + M)A(r) = 



(83) 



notice that in this way we have fixed the values of the functions Q(r,t) = e~' lEt A(r), and 
G = B(r)/A(r) appearing during the process of separation of variables. Substituting into 
flS2|) and (|8~3"D the form of the scalar and the Coulomb potentials we get 
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4- + - ) A(r) - (E - M + -(a + a') ) B(r) = 
dr r \ r 



(^-7) B(r) + (^ + M-i(a-a'))A(r) = 
Introducing the notation 

A = E + M, B = M - E, a = a + a', (3 = a' - a 
and the new variable p, related to r as follows 



P 



Dr = VM 2 - E 2 r = \f~ABi 



we have that the system of equations (p4|)-(p5|) reduces to 

Id k\ . , . ( —B a\ . 
U + P ) k[r) -{-D + -p BW = 



We shall look for solutions of the system (|S8|)-(|89l) in the form of power 



s+u a h 



A(p) = e-' 2 p 



i/=0 



u=0 



Substituting into we find 

(s + k)a — ab = 

(s — k)b — j3a = 

and 



\{s + v) + k] a v - ab u + j^K-i - &v-\ = 
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A 

[(s + v) - k] b u - j3a u + — a v -\ - b v -\ = 



D 



From (|9|)-(|3D it follows that 



s = \Jk 2 — a 2 + 



a 



a 



(95) 



(96) 



where we have dropped out ( |35D the negative root because we are looking for wavefunctions 
regular at the origin of coordinates. From (|9~3|)-(|9"5l) we have that 



[(s + v) + k] J^-P}". 



'A 

—a - [(s + v)-k]} b v 



(97) 



The series (|90| ) and (0) will have a good behavior at infinity if they terminate for a finite 
value N. Putting a^ + i = b^ + i = in (j94 T )- (|9~5|) with ^ and &at ^ 0, we arrive at 

a N D 

Substituting into (|97D, and taking into account (p6|) we get 



(98) 



[s + N)y/M 2 -E 2 = Ea + oc'M 



(99) 



where s is given by 

Finally, we obtain the energy spectrum: 



E = M 



aa 



± 



aa' 



\ a' 2 -(s + N) 2 

s + N) 2 + a 2 " \\{s + N) 2 + a 2 ) [(s + N) 2 + a 2 ] 



(100) 



Here two particular cases could be considered: a) a' = which corresponds to the Coulomb 



potential. In this case the energy spectrum reduces to: 



E = M 



1 + 



a 



-1/2 



VF 2 



a' 



(101) 



s + N) 2 _ 

where the negative root has been dropped out because it is not compatible with the relation 
(39) A second possibility is given by b) a = 0, which is the scalar V/(r) = —a'/r potential. 



In this case the energy spectrum takes the form 



E = ±M 



a 



r> 



1/2 



Vk 2 + a 



12 



(102) 



(s + N) 2 _ 

notice that in the present case states with negative energy are possible, here we do not have 
critical behavior like in the Coulomb case. 
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V. DISCUSSION OF THE RESULTS 



In this section we are going to discuss the influence of the Aharonov-Bohm potential 
and the Dirac magnetic monopole charge on the energy spectrum. Here we have mention 
that the Aharonov-Bohm as well as the magnetic monopole contributions are present in the 
expression ( |1UU ) via the factor s given by (P6"|), since the explicit form of k depends on the 



relation among m, F and g. In fact, we have that when the inequalities (p5|) or (|72D are valid, 
the expression for s takes the form 

s = yW \g |) 2 - g 2 - a 2 + a' 2 (103) 

and no contribution of the Aharonov Bohm potential is observed in ( |100| ). The values of m 
for which (|103|) takes place are given by the expression (|74]) 



It is worth mentioning that when the magnetic monopole contribution is absent, the 
inequalities given by the expression ([74]) never take place, and consequently, the expression 
for s given by eq. (|103|) is not applicable. In this case the energy spectrum can be computed 
by substituting the expression 



n+ | m - F | +1/2) 2 - a 2 + a' 2 (104) 

into ( |100|) . Analogously, the energy spectrum when F ^ and g^ can be obtained after 
substituting into (|100|) the following value of s 



n + m- F - g + 1/2) 2 - g 2 - a 2 + a' 2 (105) 

when m — F > and m — f — 2g > 0. Otherwise, when m — F < and m — F — 2g < 0, 
the value of s to substitute into (11001) reads 



n-m + F + g + l/2) 2 -g 2 -a 2 + a' 2 (106) 

Perhaps the most interesting and puzzling result of the present paper is the non dependence 
of the energy spectrum on the Aharanov-Bohm potential for a range of values given by the 
inequality ([741). Despite this phenomenon was already pointed by Hoang et afl there are 
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some discrepancies between their results and those ones present in this paper. Basically the 
problems lie on the criteria for establishing the boundary conditions and the normalizability 
of the wave functions. It is worth mentioning that since the Aharonov-Bohm contribution 
F can take non integer values, then the parameters a and (3 in the Jacobi Polynomials 
P^*'^(x) can be negative provided that a > —1, and (3 > —1. Obviously if the Aharonov 
Bohm potential is absent, from the results presented in Sec. 3 we have that a > and 
(3 > (or in the case when F is an integer). Hoang et al consider that a and (3 are always 
positive restricting in the way the range of validity of the solutions. A second point to 
remark is that we not only impose the normalizability of the wave functions but also the 
existence of the expectation value of the angular momentum operator, which is equivalent 
to say that / $t K QcLfidip < oo, in this way we have to impose that the spinor components 
£iand £ 2 presented in Sec. 3 should satisfy / ^1,2^^1,2 < °o • Regarding the assertion made 
in0 about the existence of quantum states forbidden for the Dirac particle in the presence 
of the Coulomb plus the Dirac monopole potentials, we have that the boundary conditions 
imposed on the wave function in the present article avoid such a anomalous behavior and 
consequently the spinor $ is well defined for any value of the parameters m, F and g. 
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